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We investigate the phase structure of the SUf{2) 5C/c(3) Nambu- 
Jona-Lasinio model as a function of the scalar diquark coupling strength. 
Above a critical coupling, the binding energy is sufficiently large to over- 
compensate the quark masses and a massless scalar diquark bound state 
emerges which leads to Bose condensation already in the vacuum. 
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1. Introduction 



Recent laboratory experiments with ultracold gases of fermionic atoms 
allow to investigate dense Fermi systems with their coupling strength tun- 
able via Feshbach resonances by applying external magnetic fields. After the 
preparation of fermionic dimers in 2003, now also their Bose-Einstein con- 
densation (BEC) ^j, i2j and superfluidity of these dimers has been observed 
[U Uj. Besides this strong coupling regime, for weak attractive interac- 
tions at low enough temperatures the condensation of bosonic correlations 
(Cooper pairs) in the continuum of unbound states occurs according to the 
Bardeen-Cooper-Schrieffer (BCS) theory. 

The BEC-BCS crossover is physically related [5] to the bound state 
dissociation or Mott- Anderson delocalization transition [Sj where the mod- 
ification of the effective coupling strength is caused by electronic screening 
and/or Pauli blocking effects. The Mott transition is a very general effect 
expected to occur in a wide variety of dense Fermi systems with bound states 
such as deuterons in nuclear matter [HI [9l UO] or diquarks in quark matter 
im [T2I \T3[ 114] . Below a critical temperature, bosonic correlations form 
a condensate and this transition appears as BEC-BCS crossover, which in 
quark matter is of particular theoretical interest due to the ultrarelativistic 
regime for massless (Goldstone) bosons [T6 t IT7]. 

A systematic treatment of these effects is possible within the path in- 
tegral formulation for finite-temperature quantum field theories. This ap- 
proach is rather general as it is relativistic and is especially suited to take 
into account the effects of spontaneous symmetry breaking. In this con- 
tribution we sketch the basics of this approach on the example of a model 
field theory of the Nambu-Jona-Lasinio type for a relativistic strongly inter- 
acting Fermi system, see [T8] for a recent review. These investigations are 
also motivated by the analogies of the strongly coupled quark-gluon plasma 
(sQGP) at Relativistic Heavy Ion Collider (RHIC) in Brookhaven [TP] with 
the experiments on BEC of atoms in traps. The further development of the 
approach may provide qualiative insights into the phases of QCD at high 
densities like the recently suggested quarkyonic phase [201 Ell 112] • Possi- 
ble evidence for a triple point related to this new phase comes from hadron 
production in heavy-ion collision experiments [23j to be further investigated 
at upcoming dedicated facilities, e.g., CBM @ FAIR Darmstadt, NICA @ 
JINR Dubna. 
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2. Formalism 

2.1. Model Langrangian and mean field approximation 

Our starting point is a NJL-type Lagrangian for three colors and two 
flavors, motivated from Fierz transformed one-gluon exchange 

C = qiid^-f^ - mo)q + Gg [ {qqf + {q^l^Tqf ] 

+ Gd Yl (9i75(^T2AAg^)(g^iC75r2AAg) , (1) 

A=2,5,7 

with C = i727o being the charge conjugation matrix, r = {ti,T2,T3) and 
Ti the Pauh matrices in flavor space and the anti-symmetric Geh-Mann 
matrices in color space. The parameter choice mg = 5 MeV, G^A^ = 2.1 
and A = 653 MeV reproduces the vacuum pion mass and decay constant. 

For this model Lagrangian we can give the partition function in its 
bosonized form 

r f f /3 ^2 I 2 |A|2 ^ 

Z = J VA*VAVaVTTexpl- J d\^^ — ^i^P^*'^"^' 

where diquark (A*, A) and meson (cr, vr) degrees of freedom appear as 
collective flelds instead of the quark ones which have been integrated out 
leading to the determinant of the quark propagator S in Nambu-Gorkov 
representation. For details of the further calculation we refer to [24j. By 
minimizing the thermodynamical potential, we obtain gap equations which 
need to be solved self-consistently. The corresponding order parameters 
then characterize the phase structure of the model. 

2.2. Gaussian fluctuations and polarization matrix 
We expand around the mean fleld Nambu-Gorkov quark propagator 

and obtain 



Smf = g~ ) 7 up to second order in the matrix S = ^ ^'^^^^ ^ 



Indet^-^ = Trln5^/^ + Tr (^Smf^ - ^Smf^Smf^) + ©(S^) . (3) 

The fluctuations of the collective fields can be decomposed according to: 
7r+/^ = i75T/* • vr, 6^^^ = i75T2A25*/, and their amplitudes can be arranged 
in the vector = {vr, a, 5,5*}. Performing the trace operations over Nambu- 
Gorkov, fiavor, color, Dirac and momentum space we introduce the elements 
of the polarization matrix n(A;o,k) (for details see 



^Tr (^Smf^Smf^^ = (t>i'nij4>j , (4) 
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with i,j = {it, a, 5, 5*} denoting the channels. Some matrix elements are 
pairwise equal, e.g., IIo-^ = Ils*a,^Sa = ^aS* and for real A even II^^ = 
n^*^* . Thus, the pions are degenerate, as expected for isospin symmetric 
matter. We explicitly include the mixing terms between the a and the 
diquarks in our investigation, which has been omitted in the literature so 
far \12\ [T^ [27] . Performing the gaussian path integral over the fluctuation 
fields results in an expression for the thermodynamical potential 

Q[T, /u) = -TlnZ = Indet Smf + Indet [6ij/{2Gi) - ^(w, k)] , (5) 

where Gi = {Gs,Gs,Gd,Gd}- The mass spectrum of quasiparticle modes 
can be found from the condition of the vanishing determinant in Eq. ([5]) at 
k = for cj = Wj(k = 0) = {m,r, "ig-, — fi, ms* + /x}. 



3. Results and Discussion 

We want to discuss first the vacuum case /i = T = 0, where diquark 
and anti-diquark are degenerate, the general discussion will be given in 
[26j . In the normal phase, for the dimensionless diquark coupling strength 
rjD = Gd/Gs in the range [T7] 

vr^ 3 m ^ , , 

<VD<^ = vh, (6) 



AGs (A VA^ + m2 + m2 In A+VA"+m^ ^ 2m -mo 

the polarization matrix is diagonal and thus Q for the diquarks gives 

^^^^''l'(2^{Ep + mD/2^ Ep-mD/2) ■ 

If r]D > vh the matrix is not diagonal anymore. Neglecting the mixing 
terms, we get two solutions for ([5]), namely 

' = ^^^Iw)^{E^ + mo/2^ E^-mo/2) 

^ = ^^^i Wf\^El + ^^)\E^ + mD/2 ^ E^- mD/2) 



where E^ = ^E'^ + A'^, Ep = y^p^ + m^. A Goldstone mode mo = 

solves the first equation, which in this case coincides with the gap equation 
for the pairing gap. The solution of the second equation gives a massive 
mode. The results are shown in the left panel of Fig. [1] as a function of ijd- 
In the right panel of Fig. [T]the phase structure of the model is shown for four 



zablocki'Acta'VS printed on December SI, 2009 



5 



v 

J , ^ ■ . ■ J , 

100 200 300 400 

H [McV] 

Fig. 1. Left panel: Vacuum diquark masses ruD and threshold energies iSthr = 
2v' + as function of the couphng strength r]£). Right panel: Phase diagram 
for different values of the coupling strength rjn = 0.9, 1.2, 1.5, 1.8. First and second 
order phase transitions are indicated by solid and dashed lines respectively. Dotted 
lines denote crossover transitions. The black dot for rjn = 0.9 indicates the critical 
endpoint for first order phase transitions. For the discussion see text. 

cases of coupling strenths: r/£) = 0.9, 1.2, 1.5, 1.8. While for rjo < 0.9 there 
is no coexistence of chiral symmetry breaking and diquark condensation, in 
the range 0.9 < rjD < one obtains Bose condensation of bound diquarks 
in such regions of coexistence. At > Vd ^ ^^^^^ more spectacular effect 
occurs: the vacuum state itself is a Bose condensate of diquarks! While 
this model description of a relativistic Fermi system at arbitrary coupling 
is surely of methodological interest in the context of experiments with Bose 
condensates of atoms in traps, its relevance for the discussion of the phase 
structure of QCD requires a careful analysis of the corresponding hadron 
spectrum. 
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